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Abstract 

f-*>) \ As a black hole in a binary spirals in gradually from large separation, energy 

(N; 

f^ | and angular momentum flow not only to infinity but also into or out of the 

<D : 

C/3 \ hole. In addition, the hole's horizon area increases slowly during this process. 

t— i ; 

In this paper, the changes in the black hole's mass, spin, and horizon area 

during inspiral are calculated for a hole in a circular binary with a companion 

00 ! 

body of possibly comparable mass. When the binary is composed of equal- 
ly ; 
^^ \ mass black holes that have spins aligned with the orbital angular momentum 

^H \ 

and are rapidly rotating (with spins 99.8 percent of their maximal values), it 
O 

is found that the fractional increase in the surface area of each hole's horizon 
5h '■ 

is one percent by the time the binary spirals down to a separation b of 6M 
>' 
k>H | (where M is the binary's total mass), and seven percent down to b = 2M. The 

flow of energy and angular momentum into the black holes' horizons changes 

the number of gravitational-wave cycles in the LIGO band by no more than a 

tenth of a cycle by the time the binary reaches b = 2M. The results obtained 

in this paper are relevant for the detection and analysis of gravitational waves 

from binary systems containing a black hole. 

I. INTRODUCTION 

Binary black holes are expected to be among the primary sources of gravitational waves 
for interferometric detectors 0. Since we do not have exact solutions of Einstein's equations 



that represent binary black holes in sufficient generality, we must study these systems per- 
turbatively and/or numerically. One regime in which the evolution of binary black holes is 
well understood is the early inspiral phase. In this phase, the holes' separation is still much 
larger than the binary's total mass, and post-Newtonian expansions can be used to analyze 
the system. Eventually radiation reaction drives the holes together and the post-Newtonian 
approximation fails. The binary's subsequent evolution must be studied numerically. 

While the flow of energy and angular momentum to infinity during inspiral has been 
calculated to high post-Newtonian order, to date the flow into or out of the black holes' 
horizons has not been computed except in the extreme-mass-ratio limit; and in that limit, 
it has been done to very high post-Newtonian order || (for numerical work, see e.g. ||). 
Absorption (or emission) of energy and angular momentum by the holes' horizons, while 
much smaller than emission to infinity, might still be important because extraction of weak 
gravitational signals from noisy detector output using matched filtering requires knowledge 
of the orbital evolution to very high accuracy, and black hole absorption/emission might 
affect the evolution at that level. Two purposes of this paper are to calculate black hole 
absorption/emission of energy and angular momentum to leading order in a circular binary 
with holes of possibly comparable mass, and to investigate whether it is relevant for detection 
and analysis of gravitational waves. 

A third purpose of this paper is to provide some information on the interface between 
the inspiral and merger phases of binary evolution. Numerical simulations of binary black 
holes typically begin computing at this interface and need initial data representing holes 
that have spiraled in from infinity, i.e. initial data tied to the inspiral phase and to the 
post-Newtonian expansions used to describe it. One approach to obtaining such initial data 
is given in ||. Since initial data of this sort are not yet being used, one needs to relate the 
masses, spins, and horizon areas of the black holes present in currently used initial data to 
the corresponding quantities when the holes were infinitely separated. For this purpose, it 
is necessary to know how these quantities change during inspiral. In this paper, I calculate 
the leading-order changes in the holes' masses, spins, and horizon areas during inspiral for 



a circular binary. 

Recently, Price and Whelan || have emphasized the role of angular momentum absorp- 
tion/emission by rapidly rotating black holes at the end of inspiral, when the holes are 
beginning to merge. Here I focus on the earlier stages of inspiral, when the holes are widely 
(or moderately) separated and their gravitational effects on each other can be described 
using black hole perturbation theory. 

The results obtained in this paper are actually valid for a black hole in a binary with any 
companion body (e.g. a neutron star) that is well separated from the hole. The formulas for 
the changes in black hole quantities presented here depend only on the companion body's 
mass and not on its internal structure. These formulas therefore remain valid when the 
companion is not a black hole. 

II. FRAMEWORK 

I follow the field-theory-in-flat-spacetime notation used in the literature on post- 
Newtonian expansions (e.g. 0,0) and denote 3- vectors by bold-face letters. A centered 
dot between 3-vectors denotes the usual inner product in flat 3-space; a hatted 3-vector 
represents the unit vector in that direction. 

Consider a black hole binary undergoing circular motion with the separation b between 
the holes much larger than their total mass M = Mi + M 2 , where Mb denotes the mass 
of the 5th hole. Define fj, = M ± M 2 /M and 77 = fi/M. Label the holes BH1 and BH2, and 
denote their spins by S# and horizon areas by Ab for B = 1,2. Let Sb = (Sb • Sb) 1 / 2 be 
the spin magnitudes, and define the parameters xb by Sb = XbM b (B = 1, 2). Throughout 
this paper I assume xb < 0.998 (B = 1, 2); this restriction is based on the analysis in ||. 

Define each black hole's horizon radius thb — Mb[1 + (1 — Xb) 1//2 L angular velocity 
&hb = Xb{^thb)~ 1 , and surface gravity Kb = (1 — XbY^^hb)^ 1 (B = 1, 2). Introduce the 
following Newtonian quantities for the binary: the orbital angular momentum L/v, the orbital 
angular velocity £l N = (M/6 3 ) 1 / 2 , and the relative velocity v = (M/b) 1 ^ 2 . By assumption, 



i;< 1. 

Since the black holes are widely separated, each hole has a surrounding region that 
satisfies the following properties: (i) it is far enough from the hole that gravity is weak 
there; (ii) it does not extend so far that the companion hole's tidal field varies appreciably 
in the region ||. We can place in this region an inertial coordinate system in which the 
hole is (momentarily) at rest. This region and its local coordinates are referred to as the 
black hole's local asymptotic rest frame (LARF) ||. Label the two regions around the holes 
LARF1 and LARF2. 

Usually mass and angular momentum are only defined globally in general relativity, using 
fields at infinity, since precise local definitions are not available. However, for a black hole 
well separated from its companion, one can define the hole's mass and angular momentum 
using fields measured in the hole's LARF; these definitions are inherently ambiguous |9|-|T2 



(For further discussion of the ambiguities, see Sec. |V]].) I refer to these definitions when 
discussing a black hole's mass and angular momentum in this paper. I calculate the rates of 
change of these quantities as measured in the LARF — that is, with respect to time t measured 
by an inertial observer in the LARF. When integrated over the duration of inspiral, these 
rates of change should give results exceeding the ambiguities in the definitions of mass and 
angular momentum, in order to be relevant to the analysis of initial data at the interface 



between inspiral and merger. This issue will be discussed further in Sec. [VI. 

I also consider slices of constant time t that begin in the LARF and extend into the 
black hole, intersecting the horizon in 2-surfaces that correspond to constant ingoing-time 
slices of a Kerr black hole's horizon. (Alternatively, one can consider slices that intersect a 
"stretched horizon" as discussed in ]13[ and references therein.) The rate of area increase 



of these 2-surfaces can be calculated using the results of Hawking and Hartle MM combined 



with black hole perturbation theory |L5],[T6|. The quantities dMs/dt and dSs/dt can then be 
obtained from dAs/dt using the first law of black hole mechanics dM = (k/87t)gL4 + VtndJ 
and the relation uodJ = vndM for black hole perturbation modes of angular frequency uo 
and azimuthal angular number m [|T!],|l7],[L| . (Here J refers to the black hole's angular 



momentum.) 

Throughout this paper, I focus on BH1 and the changes in its parameters. The corre- 
sponding formulas for BH2 are simply obtained by exchanging the subscripts 1 *-* 2 in the 
final results [e.g. Eqs. (PI)]. In Sec. [TTJ, I consider the special situation in which BH2 is held 
stationary with respect to BH1. The results from this artificial scenario are used in Sec. [TV] 
to analyze a circular binary with black hole spins aligned or anti-aligned with Ljy. The more 
general case of spins not fully aligned or anti-aligned with Ljv is treated in Sec. |V|. 

III. STATIONARY COMPANION 

In this section, I calculate the tidal distortion BH1 suffers when BH2 is held stationary. 
This involves solving for the Weyl tensor component ipo, which contains complete information 
about the gravitational perturbation on BH1, using the Teukolsky formalism |15|]. With ijj in 
hand, the rates of change of BH1 parameters can be calculated using the results of Hawking 



and Hartle |T^] and Teukolsky and Press [[RJ 



The first step in this process is to calculate BH2's tidal field as seen in LARF1. I 
will consider only the lowest-order Newtonian tidal field, which is approximately constant 
throughout LARF1. To calculate this field and its effect on BH1, consider first a fictitious 
Euclidean 3-space containing a single stationary body of mass M 2 at coordinate location 
(b,8o,<f)o) in a spherical polar coordinate system. The Newtonian gravitational potential at 
the field point (r, 6, </>) is given in these coordinates by 

$(r, 6, 0) = -4vr-^ £ £ (21 + l)" 1 (-) Tfafo, <f> )Y lm (6, 0) (1) 

£=0 m=—l 

for r < b. 

We are interested in the gravitational field only in a small neighborhood of the origin 
satisfying r ^ b. In particular, we would like to evaluate the body's tidal field at the 
origin, so only the I = 2 part $ ( ^ 2 - ) of $ is relevant. The (electric-type) tidal field is given by 

(2) 

Sij = $^ ■ in Cartesian coordinates. Taking these derivatives in spherical coordinates and 



evaluating in the usual spherical orthonormal basis yields the tidal field components S§§, 
£g7, Si 7 near the origin r = 0. The particular combination of relevance to us (see below) is 
in this way determined to be 

£ U - £ §§ - 2i£ §i = ^§^ E 2 Y 2m (0, 4>)Y* m (e , O ). (2) 

00 m=-2 



Here the functions 2^2m(^ 5 0) are spin-weighted spherical harmonics ITS 



Return now to the black hole binary. The region near BH1, including LARF1, can 
be described as a perturbed Kerr black hole, and so can be covered by Boyer-Lindquist 
coordinates (t,r,6,<f)). We would like to solve the Teukolsky equation |15[ in this region 



for the Weyl tensor component tpo(r,8,<j)). If we were considering a single perturbed Kerr 
black hole as the entire spacetime, the asymptotic form of ipo as rjM\ — > oo would be 



the combination £11 — £$ — 2i£§i of the external tidal field ||13f| , since ipo vanishes for 
an unperturbed black hole. In our binary system, ipo acquires this asymptotic form for 
Mi <C r <C b, i.e. in LARF1, with the tidal field £ij being that of BH2. To lowest order, 
this tidal field is exactly the Newtonian field of a body of mass M 2 at separation b, which 
was calculated above; in particular, the angular dependence of ipo in LARF1 is given by 
Eq. d|), but with 9 and (f> now representing Boyer-Lindquist coordinates, and #o and 0o 
now representing BH2's angular coordinates as seen in LARF1. Therefore, to solve for the 
perturbation ipQ on BH1, we impose the LARF1 boundary condition 

^o - ^^ E 2Y 2m (6, 01£,(*o, fa) (3) 

for Mi < r < b. 

It now remains to solve the Teukolsky equation for ?/>o with the boundary condition @ 



and an appropriate no-outgoing- wave boundary condition at the black hole horizon [ 15 1 . We 
express ipo as a sum of modes 



2 



^o= E 2Y 2m (8A)Rm(r) (4) 

m=-2 

and solve the radial Teukolsky equation for -R m (r) subject to the no-outgoing- wave boundary 
condition at the horizon. This yields the radial functions (Eq. (5.7) in Ch. VI of fT9H ) 



R m (r) = C m x^- 2 (l + x)-^- 2 F(-4, 1, -1 + 2 7m , -x) (5) 

for m/0. Here 

imxi _ r-r m 

7m ~2(l- X ?) 1/2 ' X ~2M 1 (l-x?) 1/2 ' lJ 

and F is a hypergeometric function. The m = mode can be treated separately; since a full 
treatment reveals that this mode does not contribute to the rates of change of black hole 
parameters, I ignore it here. The constants C m are determined by imposing the LARF1 
boundary condition (||); we obtain 

C m = ^^lm(lm + 1)(4t& - 1)Y*J9 , O ). (7) 

The leading-order tidal distortion of BH1 due to the presence of a stationary companion 
of mass M 2 has now been determined. This information allows us to calculate the rates 
of change of BH1 quantities using the results of Hawking and Hartle |14| . In fact, given 



the modal decomposition @, we can easily obtain the relevant rates using explicit formulas 
provided by Teukolsky and Press fl6fl . The results are dMi/dt = and 



1 ~ixi + ^xi^ 2 0o) 



dA 1 647rMfM|xJsin 2 6 (^ 3 2| 15 

~dT ~ 56 6 (l-xi) 1/2 

dSi (1 - Xi) 1/2 dA x SMfM 2 . / 3 15 2 . 

•Xi sm °o 1 - 7X1 + -rXi sm ^0 



dt 8ttxi dt 56 6 ^' u V 4^ x 4 



(8) 



Here 9 is BH2's ^-coordinate — that is, its polar angle with respect to Si — as measured in 
LARF1. 

Since the effects of only the leading-order tidal field were taken into account above, the 
expressions (§) are actually valid for any companion body of mass M 2 , not just a black hole. 
The rates (|8|) of area increase and spin-down have already been derived by Teukolsky (19 



in the extreme-mass-ratio limit, i.e. for M 2 <C Mi. The derivation I have presented above 
establishes the validity of the expressions (§) for comparable-mass black holes as well. Hartle 
and collaborators [|T^,^,^T|,|T^] have shown that the spin-down of a black hole by an external 



tidal field is analogous to the Newtonian tidal friction process in a planet-moon system. 



The results @ will be used in the next sections to obtain the corresponding formulas 
for a binary undergoing circular motion. 

IV. EQUATORIAL ORBITS 

In this section I study special configurations of the binary in which the black holes are 
in a circular orbit and have their spins aligned or anti-aligned with the orbital angular 
momentum Ljv- In these scenarios there is no precession of the angular momenta: the spins 
remain aligned or anti-aligned with L^r. As a result, the companion to each of the holes 
orbits in the hole's equatorial plane; more precisely, the external tidal field seen by each of 
the holes rigidly rotates about an axis parallel or antiparallel to the hole's spin axis. 

In Boyer-Lindquist coordinates (t,r,9,(fi) centered on BH1, with Si along 9 = 0, the 
t- and ^-dependence of the companion's tidal field enter in the combination <ft — Qt. The 
rotation rate Q of the tidal field as seen in LARF1 is to leading order Q = (L^ ■ Si)f2jv, 
where Ljy • Si = +1 (—1) for a prograde (retrograde) orbit. The first correction to this 
expression for Q is 0(v 2 ) higher (see Eq. (3.12) in [|J), and will be ignored in this paper. 

A. Instantaneous rates 

In the rigid 0-rotation case, simple formulas given in Eqs. (7.21) of [|13j (and reproduced 
below) specify the rates of change of black hole quantities in terms of a horizon integral / 
that depends on the particular perturbing gravitational fields present: 

dSl (n-n*i)i, ^ = n^ = w-n m )i, 



dt " dt dt 



(9) 



In terms of ingoing Kerr coordinates (V, r, 9, (ft) (see e.g. |22| for a definition), / is an integral 
of a function of 9 and (0 — QV) over a constant-^ slice of the horizon. Since 0-rotations are 
isometries of the horizon metric, I is independent of V. 

8 



Consider an expansion of I in powers of M±Q, which is 0(v 3 ) and hence much smaller 
than 1. The zeroth-order part I = I\ n=0 is independent of Q and, in our situation of 
binary black holes, can be easily obtained from the results for a stationary companion. 
From Eqs. (|9|), we have St = —Qnih, where an overdot indicates a time derivative. 

But Q = corresponds to a stationary companion, and in this case we have an explicit 
expression for Si in Eqs. fl8|). Equating 5i in Eqs. (f|) to — QhiIo yields 

Wo) = l ~^MlMl sin 2 6 (l - \ X \ + ^xi sin 2 O ) , (10) 

where 9$ = rr/2 for the equatorial orbits considered in this section. The general expres- 
sion (0) with a wider range of values for 9q will be used for non-equatorial orbits in the 
next section. Since the first correction to Iq in the expansion of I in powers of Mitt is 
O(MxQ) = 0(v 3 ), I will approximate I by Iq throughout this paper. 

Note that Eqs. (§) are, strictly speaking, valid only for constant rotation rates Q. In our 
situation, radiation reaction drives the binary together and so Q changes during inspiral. 
However, the timescale for these changes is the inspiral timescale r- ms ~ bv ~ 6 , where "~" 
means "is of the order of"; this is to be compared to the timescale k{ 1 on which the 
divergence and shear of the null generators of the horizon probe the future |]il||iT| , |i"5[ . By 



assumption, xi is l ess than or equal to 0.998; this implies k{ 1 < 34Mi, so k{ 1 is much smaller 
than r ins . Therefore Eqs. (|9|) are valid in our binary system to a very good approximation. 
The various timescales of interest to us will be discussed in more detail below. 

Note also that Eqs. (||) [and, in addition, Eqs. (|TTD, (^l|), and (|22|) below] are valid only 
when integrated over time intervals much longer than n^ 1 (see the discussion in Sec. VI.C.ll 
of |13|| ). In this paper, I am interested in integrating these equations over the entire inspiral — 
that is, over time intervals of order r ins — so this condition is certainly satisfied. 

After putting I (ti/2) and Q = (Ljy • Si)f2jv into Eqs. (|9|), we obtain the following rates 
of change of BH1 quantities for a circular orbit with spins aligned or anti-aligned with L^: 

^ = (fi-fim)/o(7r/2) 



'dJ\ v 5 /M M ., ., 



-Xi 



1 + 2(L w -S 1 )[l + (l-x?) 1/2 



eft 






fi(fi - n Hl )I (ir/2 
'dE\ - 5 ' "' ^ ' 

/ JV 
87TK^ 1 (fi-fi m ) 2 / (7r/2) 

64vrM 1 5 M 2 2 (l + 3x?) 






Mj 3 

IT f . 



Mi 



56 6 (l-xi) 1/2 



Xi-2(L JV -Si) 



1 + fl 



xl) 1/2 



Mi 3 
— tr 
M 



(11) 



In these formulas, the Newtonian quadrupole expressions for energy and angular momentum 
flow to infinity are [f23pl 



'dE 
dt 



32 



V 2 v 10 , 



dJ 

dt 



32 



rfhlv' 



(12) 



N V ""' J N 

where v = (M/b) 1 ^ 2 and i] = MiM 2 /M 2 . Note that energy and angular momentum absorp- 
tion/emission by a rotating (non-rotating) black hole is 2.5 (4) post-Newtonian orders below 
the quadrupole emission (|i~2"l) to infinity, as first derived in the extreme-mass-ratio limit by 
Poisson and Sasaki |25j and Tagoshi, Mano, and Takasugi ||. The rates of change for BH2 
are obtained by exchanging the subscripts 1 <-> 2 in the formulas (ll]). 



The energy absorption/emission rate Mi given above agrees in the limit M 2 /M — > with 
the lowest-order expression obtained by Tagoshi, Mano, and Takasugi 0. Those authors 
have calculated this rate in the extreme- mass-ratio limit, for a circular equatorial orbit, to 
much higher order in v than I have done here. However, their results are not applicable to 
comparable-mass binaries, while the formulas ( |Tl~D are. 

The expressions (|TT| ) are valid even if BHl's companion is not a black hole, provided the 
companion's mass is substituted for M 2 . 



B. Total changes during inspiral 

In this subsection, I integrate Eqs. flTT| ) to calculate the total changes in M±, S\, and A\ 
during inspiral. I take into account only the leading-order Newtonian effects of radiation 



10 



reaction when computing orbital decay; given this approximation, the orbital separation b 
evolves as p4] , |2"2| 



b(t) = 6 (1 - i/T"o) 1/4 , 



(13) 



where r = (5/256)6g(;uM 2 ) _1 . I also ignore all post-Newtonian corrections to the orbital 
angular velocity Qn- 

It is convenient to parametrize the orbit by separation b instead of time t. The total 
change in a parameter, say S±, from infinite separation to separation b is denoted ASi(b) 
and is calculated by integrating Eqs. fllTl) . As a first approximation, the quantities Mb 
and Sb (B = 1, 2) on the right-hand sides of Eqs. fllTf) can be considered constants during 
inspiral. The reason is that the timescales for evolution of Mb and Sb are much longer 
than the inspiral timescale r ins ~ bv~ e . Indeed, the timescale for evolution of the masses is 
tm ~ Mb/ Mb ~ bv~ 13 , and for the spins is t$ ~ Sb/ 'Sb ~ bv ~ 10 . So tm 3> Ts ^ r ms and 
we can safely treat M B and Sb (-8 = 1,2) as constants on the right-hand sides of Eqs. ( |TT|) 
when integrating over inspiral. 

With these approximations, the normalized changes in BH1 parameters from infinite 
separation to separation b are 



AS 1 vMx 2 

-^-( fo )^TTr( 1+: >v,) 
AMi 

AAi 



Ai 



(6) 

(6) 



4M 

4 U/ V 
»jA/?(l + 3*f) 



Xi 

' 4 



/ /lf\ 

Kt) +(^-Si)[i + (i-x?) 



2x1/2 



2Mi /M\ 7/2 



— 



7M \b J 



+ 3X?) 



:Lv Sl) f( T ) + 



i+d 



x 2 ) 1/2 



Mi /My 



5M \ b J 



8Mr 



Hl\ 



+ 



4 / rm 

5 VM 



Xi) 1/2 



M 



6 / 



,f A x8xir m /M\ 7/2 



(14) 



where r#i = Mi[l+(1— Xi) 1 ^ 2 ]- To evaluate these changes, one can put into the formulas fll4|) 
the values of Mi, Si, and Ai at infinite separation or, for that matter, at any separation 
much larger than M, because the changes in these quantities during inspiral are small. 
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Xi b/M = 100 b/M = 20 b/M = 6 b/M = 2 

9. x lO" 10 2. x 10~ 7 2. x 10~ 5 

0.5 -7. x 10~ 7 -2. x 10~ 5 -2. x 10~ 4 

0.998 -3. x 10~ 6 -8. x 10~ 5 -8. x 10" 4 -6. x 10~ 3 

TABLE I. Normalized change ASi/Mf in spin evaluated at 6/M=100, 20, and 6 for an 
equal-mass binary with Ljv ■ Si = 1. For rapidly rotating holes (xi = X2 = 0.998), this change is 
also evaluated at b/M = 2. 

Once again, the changes for BH2 are obtained by exchanging the subscripts 1 <-> 2 in the 
expressions (|14|). 



The normalized parameter changes (H), evaluated at different stages during inspiral, 
are displayed in Tables |-fT| for an equal-mass binary (Mi = M 2 ) with L^v • Si = 1. Since 
a binary composed of slowly rotating black holes is expected to be undergoing a transition 
from inspiral to merger by the time it reaches b = 6M, the endpoint of integration is chosen 
to be b = 6M when xi — and 0.5. For rapidly rotating holes (xi = X2 = 0.998), the 
endpoint is chosen to be b = 2M. The assumption M <^C b is not valid at and near these 
endpoints. The results presented here are most accurate in the early stages of inspiral, when 
the black holes are widely separated, and are a rough estimate of the true parameter changes 
in the late stages of inspiral. 

C. Effect on orbital evolution 

The orbital evolution of binary black holes is affected by the absorption/emission of 
energy and angular momentum by the holes. In particular, the number of orbits — and 
hence the number of gravitational- wave cycles emitted to infinity — changes when black hole 
absorption/emission is accounted for. To estimate this effect, let us consider a circular, nearly 
Newtonian binary, with spins aligned or anti- aligned with L^r, that is losing orbital energy 

12 



Xi 


0.5 
0.998 



b/M = 100 

3. x 10~ 13 

-2. x 10~ 10 

-9. x 10~ 10 



b/M = 20 

1. x 10~ 9 

-5. x 10~ 8 

-2. x 10~ 7 



b/M = 6 

4. x 10~ 7 

-3. x 10~ 6 

-2. x 10~ 5 



b/M 



-6. x 10~ 4 



TABLE II. Normalized change AM\/M\ in mass evaluated at 6/M=100, 20, and 6 for an 
equal-mass binary with L^ • Si = 1. For rapidly rotating holes (%i = \2 = 0.998), this change is 
also evaluated at b/M = 2. 



Xi 



b/M = 100 



b/M = 20 



b/M = 6 



6/M 





0.5 
0.998 



6. x 10 -13 
2. x 10~ 7 
5. x 10~ 5 



2. x 10" 9 
5. x 10~ 6 
1. x 10~ 3 



8. x 10~ 7 
4. x 10~ 5 
1. x 10~ 2 



7. x 10" 2 



TABLE III. Normalized change A.A1/A1 in horizon area evaluated at 6/M=100, 20, and 6 for 
an equal-mass binary with L^r ■ Si = 1. For rapidly rotating holes (xi = X2 = 0.998), this change 
is also evaluated at b/M = 2. 
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and angular momentum to infinity via Newtonian quadrupole radiation (0), and to the black 
holes via tidal interaction as specified by Eqs. fll"T|). Since M B = QSb = (Ljv ■ Sb)&nSb 
(B = 1,2), circular, nearly Newtonian orbits remain circular. Therefore the evolution of 
the separation b(t) is determined by setting the rate of change of Newtonian orbital energy 
(given by E orh = —M 1 M 2 /2b) to the rate of energy loss to infinity and to the holes: 

dE orh M 1 M 2 db (dE\ dM 1 dM 2 



dt 2b 2 dt \dt J N dt dt 



(15) 



where (dE/dt)^ is given in Eqs. (|12D and Mi in Eqs. (|Tl|) [with M 2 obtained by exchanging 
the subscripts 1 <-► 2 in Eqs. (|TT|)1. 

The number of gravitational- wave cycles N± emitted to infinity from initial time t{ to 
final time tf (corresponding to separations foj and bf) is given by 

/•*/ , n N 1 rh dt (M\ l/2 



where dt/db is determined from Eq. (|T5l). This number is to be compared with the number 
of cycles N 2 obtained by ignoring black hole absorption/emission of energy and angular 
momentum, i.e. by setting E orh equal to —(dE/dt) N . The difference AiV = Ni — N 2 measures 
the effect of black hole absorption/emission on the binary's orbital evolution. 

The values of AiV obtained by setting bi to be the separation at which the gravitational- 
wave frequency is 10 Hz (the low-frequency end of the LIGO band), Xi an d X2 to be 0.998, 
and the spins to be aligned with h N (i.e. Ljv • Si = h N ■ S 2 = 1) are displayed in Table [TV] 
for various choices of total mass M (in units of a solar mass M & ) and mass ratio M\jM 2 . 
In the table, the numbers without parentheses are obtained by setting bf = 6M , and those 
with parentheses by setting bf to be the larger of 2M or the separation at which the wave 
frequency is 1000 Hz (the high-frequency end of the LIGO band). For non-rotating black 
holes (xi = X2 = 0), the corresponding values of AiV (with bf = 6M) are all less than 10~ 2 . 

The values of AiV in Table [IV] indicate that black hole absorption/emission of energy 
and angular momentum during inspiral may not be an important effect for the detection 
(by LIGO and VIRGO) and analysis of gravitational waves from comparable-mass black 
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M{M G ) M x /M 2 = 1 M 1 /M 2 = 2 M 1 /M 2 = 4 

5 0.07 (0.07) 0.11 (0.11) 0.23 (0.24) 

20 0.05 (0.07) 0.07 (0.10) 0.16 (0.22) 

50 0.03 (0.06) 0.05 (0.08) 0.11 (0.18) 

TABLE IV. Change AN in the number of gravitational-wave cycles due to black hole ab- 
sorption/emission, for various values of total mass M and mass ratio M\/M 2 . The initial sep- 
aration is such that the wave frequency is 10 Hz and the spins satisfy Xi = X2 = 0.998 and 
J-iN • Si = Lat • S2 = 1. The numbers without parentheses are for a final separation bf of 6M; 
those with parentheses are for bf equal to the larger of 2M or the separation at which the wave 
frequency is 1000 Hz. 

holes. Indeed, post-Newtonian corrections to the equations of motion and energy loss have 
far greater influence on the number of wave cycles emitted by the binary f||7]]. It should be 
noted, however, that black hole absorption/emission could have a much larger impact on the 
orbital evolution of rapidly rotating holes when they are beginning to merge, as suggested 
by Price and Whelan ||. They have presented models in which the tidal torque that results 
from black hole absorption/emission of angular momentum plays a crucial role in the late 
stages of binary evolution (see Fig. 1 in ||). The perturbative methods used in this paper 
(based on wide separation of the binary) are not valid in the close limit analyzed in || . 

It has been pointed out by Hughes |J that in the extreme-mass-ratio limit, black hole 
absorption/emission can strongly influence the binary's orbital evolution and is an important 
effect for LISA. 

V. NON-EQUATORIAL ORBITS 

In general, binary black holes are not expected to have spins aligned with the orbital 
angular momentum L^r- This misalignment causes the spins and orbit to precess in a com- 
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plicated way due to spin-orbit and spin-spin coupling 0,0] . Each black hole's companion 
is not in general confined to the hole's equatorial plane, and so the formulas in the previous 
section are not applicable. However, for orbits suitably close to the equatorial plane (see 
below for details), one can imagine using an approximation scheme in which at each instant 
the companion's #- velocity is ignored; that is, the companion is taken to be rigidly rotating 
in the ^-direction at each point on the orbit. The changes in black hole parameters can then 
be calculated by putting the instantaneous 0- velocity into the rigid 0-rotation formulas © 
at each point on the orbit. In this section, I construct such an approximation scheme. 

A. Description of orbit 

The evolutions of the spins and orbit are described by the equations [p6,7H 



S B = ftgi x S B , t N = O orb x L^ - —ri 2 Mv 7 i N , (17) 



for B = 1,2. The orders of magnitude of the precession frequencies are fi spin ~ t> 3 6 _1 and 
^orb ~ v 4 b~ l . Since the Newtonian angular velocity is Q^ ~ vb~ x , both f2 spin and f2 or b are 
much smaller than Vt^. This means that, over a few orbital periods, Ltv(£) and Sb(£) do 
not change much due to precession. Thus, the companion's orbit as seen in LARF1 is to 
a good approximation confined to a single plane with normal vector n = Ljv(t) along the 
instantaneous direction of the orbital angular momentum, on timescales of a few orbital 
periods. 

In this subsection, I analyze the trajectory of a particle in a planar, circular orbit of 
arbitrary orientation in a fictitious Euclidean 3-space in terms of spherical coordinates. 
This information will be used to specify the rotation rate and orientation of the companion's 
tidal field as seen in LARF1. Denote the particle's radial coordinate by b, its constant (non- 
negative) angular velocity by u, and the normal to its orbital plane by n. The angle of 
inclination of the normal with respect to the z-axis is denoted 8 n ; so cos8 n = n • e^ = n z . 
Assume the orbit is centered on the origin, so the particle's position X(t) at time t is given 
by a rotation -R(n, uot) about n, by an angle out, of the initial position Xo. 
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In Cartesian coordinates, the particle's trajectory is given by X(t) = X cosc<j£ + 
(n x Xq) smut. In terms of the particle's angular coordinates 9{t) and 4>(t), X(t) is equal 
to b[sin9(t) cos 4>(t), sin 9(t) sin<fi(t), cos 9(t)]. I choose the initial position to be in the equa- 
torial plane, i.e. Zq = Xo ■ e^ = 0. This choice does not affect the orbit-averaged quantities 
I calculate later in this section. 

The angular functions 9(t) and <p(t) can now be expressed in terms of n and u using the 
relations above. The quantities of interest are sin 2 9(t) and <fi(t), which are determined to be 

sin 2 9{t) = 1 - sin 2 9 n sin 2 out, 0(t) = UC ? S ^ . (18) 

sin 9[t) 

B. Approximation scheme 

Return now to our black hole binary, and go to Boyer-Lindquist coordinates (t, r, 9, </>) 
in LARF1. The companion's trajectory as seen in LARF1 will be described (to lowest order 
in v) by angular functions 9{t) and (f)(t) given by the expressions ( |T8"D with u replaced by 
Qn and 9 n now referring to the angle of inclination of Ljv(i) with respect to Si(t), that is, 
cos9 n = Ljy(t) • Si(t). After these substitutions, we have 

sin 2 9(t) = 1 - (l - \t N (t) • S^)! 2 ) sin 2 ^, 0(t) = ggMfl^i® , (19) 

V L J / sin 9(t) 

Since 9 n is now time dependent, these expressions are meaningful only when used to calculate 

orbit-averaged quantities. 

Consider the regime in which sin 2 9{t) and <fi(t) are slowly varying; more precisely, require 

them to be approximately constant on the timescale k^ 1 associated with the horizon. As 

noted before, this is the timescale on which the null generators of the horizon probe the future 



I^ , ^7|JT3| . The teleological behavior of the horizon is, however, exponentially limited; that 
is, the influence of future events on the horizon decays exponentially in time, with decay 
rate K\ (see, e.g., the discussion of teleological Green functions in fHifl). We thus require 



and 



-i 



1/sin 9)d{sin 9)/dt [which are the same to leading order by Eqs. flT9|)] to be 



only several times larger than n 1 , rather than orders of magnitude larger. 
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By assumption, Xi < 0.998, so k 1 1 is less than 34Mi. Our requirement can then be 
expressed as 

34aM^ < , 2 ff .^1 „ (20) 

sin U n \sm2\l]sft\ 

for all t, where a is a number roughly in the range 2-4. A sufficient condition for this 
constraint to be satisfied is cot 2 9 n > "iAaMiVt^- This requires Lat(£) to be near one of the 
polar axes ±Si(t), which correspond to 9 = 0, n; or, equivalently, the orbital plane must be 
near the equatorial plane. 

We are interested in separations as small as b = 6M, so VIn can be as large as (6 3 / 2 M)~ 1 . 
For this reason, I impose the constraint cot 2 9 n > 34a6 _3//2 and set a to be approximately 
3, obtaining the approximate constraints < 9 n < 7r/9 or 87r/9 ^ 9 n < n. In other words, 
Ljv(t) is within 20-degree cones around the polar axes, or, equivalently, the inclination angle 
of the orbit with respect to the equatorial plane is less than or (approx.) equal to 20 degrees. 

For the approximation scheme in this section to be valid, we require further that in the 
horizon's reference frame, the external tidal field should rotate primarily in the 0-direction 



< 



n 



in 



The 



and not significantly in the ^-direction. More precisely, we require 
rates of change presented in Eqs. ( |2~T|) and (|22"D below are subject to this condition. For most 
values of Xi, this condition is automatically satisfied throughout inspiral (down to b = 6M). 
Even if it is not satisfied at some point during inspiral, the restriction on 9 n discussed above 
ensures that the effect of the ^-motion, when integrated over inspiral, is negligible compared 
to that of the 0-motion, for almost all values of Xi- 

With the above restriction on 9 n , we can at each instant take sin 2 9{t) and 4>{t) to be 
constant relative to the horizon timescale /% x , and apply the rigid 0- rotation formulas @ 
with the instantaneous values 9(t) and <fi(t) put in. This yields 

dSl r 0(t) - Q H1 ] I o [0(t)], ^ = i(t) [<j>(t) - Q H1 ] I [9(t)], 



dt 

%TT~dT 



<j>{t)-n H1 ] i [9(t)], 

(21) 



where 6{t) and <p(t) are given by Eqs. ([OJ) and Iq by Eq. fllPf) 



C. Orbit-averaged quantities 



Next I would like to average these rates of change over an orbit assuming the binary's 
masses, spins, separation, and orbital angular momentum are approximately constant over 
an orbital period. This assumption is justified by the following ordering of the relevant 



timescales: Q^ <C Qs P - m , fi OTb <C r ins <C ts <C Tm- We can therefore take all the quantities 
on the right-hand sides of Eqs. (|2~lf) except 9(t) and <fi(t) to be constant, to a good approx- 
imation, when averaging over an orbit. Denote orbit averages by angular brackets (); so, 
for example, (Si) = (Vt N /2-n) Jq n Sidt. Plugging the expressions (fL9|) into Eqs. (|2l|) and 
performing the orbit averages (as defined above) yields 

/dSi\ 



rm 



dt 



' dMi 
dt 



m6 MlMi(l6(l + 3 X i){2Q N Afi(t) - Q m [N*(t) + l]} 

+ i5x? [M 2 (t) - 1] [m N Mi{t) - n m [3A/?(t) + 1] } 

2 -^MlM 2 2 n N Ni(t) (2n N (A - 3 Xl )sign[M(t)] - 8ft H1 (l + 3*?) 
+ 15 X l{2n N Afi(t) + n m [l - Af?(t)] } 



'dAA 87ir 2 H1 MfM$ 



dt 



16 - 12x?){^/i [A^ 2 (t) + l] - An m n N Afi(t) 



~ hb^l-xlf' 2 
+ 20^ \Mi(t)\) + ISxlfam [3A^ 4 (t) + 2A/" 2 (£) + 3 

- 8fi m ^M(t) [AT 2 (t) + l] + 8fi^/V?(i)}Y 

(22) 

where A/"b(£) = Ljy(t) • S#(£) for B = 1,2. The corresponding expressions for BH2 can be 
obtained by exchanging the subscripts 1 <->• 2 in Eqs. (|22]) . Note that these equations are 
valid only for A/s(£) suitably close to ±1, as discussed above. The formulas (^) can be 
applied to a black hole in a binary with any companion body (e.g. a neutron star) that has 
mass M.2 and is well separated from the hole. 
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Numerical integration of Eqs. (^) using the 2.5 post-Newtonian equations of motion for 
spinning bodies ( J7| and references therein) yields results comparable to those in Tables §-fT|. 

VI. DISCUSSION 

Having obtained the leading-order changes in a black hole's mass and spin during inspiral 
[see Eqs. (|14"D1, we must check whether these changes exceed the ambiguities inherent in the 
definitions of mass and spin ||. Denote by 8M and 8S the magnitudes of the mass and spin 
ambiguities. From Eqs. (1.8) in ||, 

SM ~ >£, 5S ~ "£ (23) 

n 2 n 2 K J 

where M and L are the mass and size of the (isolated) body in question, and 1Z is the 
external universe's radius of curvature. For a black hole in a binary, say BH1, L ~ M x and 
TZ 2 ~ b 3 /M 2 . This implies 

5M 1 8 Si Mi /M\ 



rxy 



Mi M x 2 ' M \b J 



V^r - . (24) 



From Eqs. (0), the changes AMi and AS'i from infinite separation to separation 6 are 
AMi /MA 2 /M\ 7/2 A5i M l (M\ 2 

~m7 ~ n m J ItJ ' mT ~ ^ m It j ■ (25) 



So, at separation 6, we have 

AMi Mi /M\ 1/2 AS 1 ! b 



— 



(26) 



(JMi M \b J ' <J5i M 

We conclude that | A^Sil exceeds the ambiguity 5 Si in the definition of spin, but |AMi| does 



not rise above 8 Mi. Note that the concept of tidal work is unambiguous |L0|-|I2|]. 

When analyzing initial data that contain a black hole and represent the interface between 
inspiral and merger, one can define and calculate the hole's mass and spin in different ways, 
giving different answers corresponding to the ambiguities 8M and 8S discussed above. Since 
8M is larger than |AM|, the hole's mass can be considered constant during inspiral to the 
same level of accuracy as used in defining mass. On the other hand, |AS| exceeds 8S, so 
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the hole's spin cannot be considered constant; however, as Table | indicates, the changes in 
spin are small during inspiral. 



The results of this work — in particular, Eqs. (|bl| ) and (|22| ) — can be used to relate the 
spin and horizon area of a black hole in a particular initial data set to the spin and horizon 
area the hole had when infinitely separated from its companion. 
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